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> \ 1 Introduction 
m 
in 

Q ■ The study of early universe cosmology demands the use of particle physics, because Hub- 
^ . ble expansion means that the universe began at a very large temperature and density. By 
! looking back to time scales on order seconds and temperatures on order IMeV, and using 
^ \ QED, nuclear physics, low energy weak interaction physics, and nonequilibrium kinetics, it 
has been possible to predict the primordial abundances of elements. It is natural to try to 
push this success back to shorter times and higher energies, where one finds a sequence of 
interesting phase transitions: the QCD phase transition, which it may soon be possible to 
r-| ! study through heavy ion collisions; the electroweak phase transition, which may be respon- 
sible for baryogenesis; and perhaps even GUT phase transitions, which may have bearing on 
k>( ; baryogenesis, inflation, and vacuum stability. 

\^ • Our attempts to understand these phase transitions, and thermal physics at these energy 



scales more generally, have been hampered for two reasons. First is an uncertainty about 
the underlying physics, in the electroweak case and especially in the GUT case. The second 
is the infrared problem of interacting, bosonic thermal field theories. Even at weak coupling, 
if one is interested in sufficiently long wavelength phenomena, mutual interactions between 
light bosonic species cause ordinary perturbation theory to break down. Perturbation theory 
is inevitably useless in the high temperature phase of Yang-Mills Higgs theories when one 
considers length scales of order l/gf^T, and for shorter length scales this problem invariably 
enters calculations at some perturbative order, because of interactions with modes at this 
scale. 

This problem was first noted by Linde |jl|, and has since been understood as arising from 
the essentially 3 dimensional nature of the thermodynamics of the infrared, bosonic fields. 
The partition function describing equal time, equilibrium thermodynamics, Z = Tre~^^, 
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can be written as a Euclidean path integral where "time" is periodic with period (3 = 1/T. 
At weak coupling, the decoupling theorem allows one to integrate out Fourier modes with 
nonzero wave number in this direction, the nonzero Matsubara frequency modes[^. This 
leaves a three dimensional theory, which is superrenormalizable, and therefore well behaved 
in the ultraviolet but potentially poorly behaved in the infrared. The couplings of this 
theory have nonzero engineering dimension, which sets a natural scale where the perturbative 
expansion becomes unreliable. The physics of interest often lies in this region. 

The solution to this problem is to perform the integration over the nonzero Matsubara 
frequencies analytically, reducing the theory to the simplest sub-theory which contains the 
correct problematic infrared behavior. Then one studies the 3-D theory nonperturbatively, 
for instance on a lattice. The lattice treatment also turns out to be particularly tractable 
because of the reduced dimension and the superrenormalizability. 

In the case of the electroweak phase transition this program has been pushed nearly to 
completion. The dimensional reduction has been performed 0, and the match between 
the 3-D continuum theory and a 3-D lattice theory has been carried out so that they will 
have the same small lattice spacing a limit This matching procedure can be exact 

because of the superrenormalizability of the 3-D theory. Several groups have studied the 
resulting system numerically 0, 0, ||, ^ |T0[. There has also been a very general analysis of 
the dimensional reduction step [|ll|] , which has been applied to the supersymmetric extension 
of the theory as well [Q. 

Numerical studies of the lattice theory display quite strong 0(a) systematic errors, which 
at first were removed by extrapolation [Q. Except for an additive correction to and the 
value of the Higgs mass parameter at the equilibrium point, all of these errors arise from one 
loop corrections to the lattice continuum match and can be absorbed into a renormalization 
of the couplings and wave functions of the lattice theory. In particular such a renormalization 
removes all 0{a) errors in the determination of observables which describe the strength of 
the phase transition. The details of this 0{a) improvement step, in the case of SU(2)xU(l) 
Higgs theory, are worked out in |jl3|. (There is a numerically unimportant combinatorial 
error there, and two in Appendix C. We will correct them here.) 

Recently there has been interest in applying the dimensional reduction program to several 
other theories. It may be possible to understand the thermodynamics of hot QCD above the 
chiral phase transition in terms of SU(3) -|- adjoint Higgs theory Superconductivity 
is described by a Landau- Ginsburg model which is 3-D U(l) -|- Higgs theory, which also 
describes the Kibble mechanism for cosmic string formation. It has been the topic of recent 
study ||16|. If the stop squark is light, then it may be necessary to study the supersymmetric 



Standard Model by dimensional reduction to SU(3) xSU(2) theory with a fundamental Higgs 
field for each group [|17|. It will also be necessary to use dimensional reduction and a 
numerical study of the reduced theory to understand the details of GUT phase transitions, 
such as the SU(5) breaking phase transition [1 



In each case it is best to perform both the dimensional reduction calculation and the 
lattice to continuum matching calculation before undertaking extensive numerical work. 
In most cases the dimensional reduction step has been performed; otherwise the generic 
rules derived in can be used. Recently, Laine and Rajantie have computed the lattice 
continuum relations for scalar masses and the order parameter for 3-D SU(A^) gauge 
theory with either a fundamental or an adjoint Higgs field, at two loops ||T9|. This is sufficient 
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to ensure the lattice model has the correct small a limit. It seems profitable also to extend 
the 0{a) matching between the lattice and continuum theories to the general case, so that 
measurables relating to the strength of the phase transition can be computed up to O(a^) 
corrections. We do this here. 

Here is an outline of the remainder of the paper. In Section |^ we will explain the general 
idea of the 0{a) improvement, and why it is possible. Section |^ will present the coefficients 
of the 0{a) corrections for SU(A^) + adjoint Higgs + fundamental Higgs, with comments 
on how to remove either Higgs field. Section ^ shows what to do when the gauge group is 
SU(A^)xSU(M), each with Higgs fields. The last section is the conclusion. There are also 
two appendicies; the first estimates O(a^) errors by computing a particularly simple and 
important subset of graphs, and the second gives numerical values in particular theories and 
step by step instructions for building 0(a) improved lattice actions. 



2 0{a) improvement: the general idea 

Considerable attention has been given to the problem of improving the match between lattice 
and continuum regulated quantum field theories. Most attention has focused on the case of 
4-D theories, which are relevant for vacuum field theory. For instance, in 4-D pure SU(A^) 
Yang-Mills theory with no fermions, the continuum action is 

j d^x^TiF^^F^,, (1) 

where F^^ = F^JT"" = {i/ g^)[D^, D^], and = - igcT^'A^^ is the the fundamental 
representation covariant derivative. Here it is understood that divergences are to be removed 
in the MS scheme at a renormalization point fi and that Qc refers to the value at that 
renormalization point. 

A corresponding lattice action is 

4l«'EETr(l-P.,(a:)), (2) 

where Pij{x) is the product of link matricies around a square in the i,j direction and starting 
at the lattice site x. 

What we want is a match of this lattice theory to the continuum theory such that they 
produce the same infrared effective theory, up to small and controlled errors. The general 
philosophy of the matching procedure was developed by Symanzik years ago [^. If the 
theory were strictly linear then the matching procedure would be trivial; setting gf = g"^, we 
could reproduce the behavior of the continuum theory up to O(a^) errors associated with 
nonrenormalizable operators. However, the nonlinearity of the theory makes the infrared 
behavior depend on the ultraviolet, so if we change the UV modes then it modifies the IR 
effective theory. Changing a is changing the UV modes, so the relation between g"^ and gf 
will be more complex. But the relation between the theories should be analytic (or at least 
asymptotic) in the coupling when it is small. So the relation can be expanded; 

9i = 9c{^ + ci{af^)9c + C2{an)gt + ...). (3) 
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To match the theories at 0(a'') we need to know the values of ALL of the coefficients q. 
To match beyond 0(a°), say at O(a^), we must write down a lattice action with enough 
parameters to separately tune the tree level values of all O(a^) nonrenormalizable opera- 
tors consistent with cubic symmetry, and the appropriate coefficients can each be written 
as an expansion in as above. In theories with fermions there are also potential 0{a) 
nonrenormalizable operator errors. 

In the early 1980's it was hoped that the series such as the one in Eq. (|^) would be rapidly 
convergent, and perturbative computation of the first few terms would give a good match 
between lattice and continuum theories In fact, the lattice perturbative calculations 



are very complicated, and the convergence of the perturbative series is not very good. It is 
more practical to make nonperturbative measurements in the lattice theory and find gf by 
matching their values to continuum (perturbative or measured) values. It is also possible to 
go beyond an 0{aP) matching; for instance, Liischer's group has recently been able to match 
lattice Yang-Mills theory with Wilson fermions to the continuum theory, using an improved 
lattice action and nonperturbative matching conditions to remove all 0{oP) and 0(a) errors 
at the nonperturbative level . 



The idea in 3-D bosonic Yang-Mills Higgs theory is similar, but because we are in a lower 
dimension, the behavior of the expansion is fundamentally different. The continuum action, 
say, when there is a fundamental complex Higgs field, is 

j d^x (^TrF,,F,, + {D4)^D4 + m^^^V + A3e(0V)') • (4) 

The dimensionalities of the fields and couplings are now [0] = energy^/^, [A\ = energy 
[g^] = [A] = energy, and [m|] = energy^. Note that the coupling constants are dimensionful. 
In terms of the 4 dimensional couplings, g^ = gfT at tree level.0 
A lattice version of this theory has as its action 

E j ^ E Tr (1 - P^J) + % E - Uiix)(t>ix + - t/,(x)0(x + t))) + 

X [ Qzfi i<j i 

+mlZ^(j)\xmx) + XsiZl I . (5) 

If the theory were noninteracting then the lattice-continuum relations would be 

gli =9lc^ Z^ = l, As/ = Age , mil = , (6) 



and the measured lattice value of the operator insertion would match the continuum 

value after lattice zero point fiuctuations were corrected for. (The only operator insertion 
we will consider in this work is (0^)-) 

The corrections to these relations must again be analytic in g\^ A3, and mg. But all of 
these quantities are dimensionful, and their dimensionality will be balanced by powers of a. 
So, for instance, 

2 2/1 I 2, 42, 2\2| \ /^7\ 

9?.i = 93c + cifi-acfl + C25'3ca + cag^c^ca +...). (7) 



^Alternately one may maintain the 4-D normalization of the fields and couplings, but then a dimensionful 
factor, 1/T, appears in front of the action, so T always appears along with or A in the loopwise expansion. 
This leads to the same parametric arguments we present here, as each appears with a T which has to be 
balanced by an a for dimensional reasons. 
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The appearance of positive powers of a in this equation is because of superrenormahzabihty. 
Because all couplings have positive engineering dimension, the theory can be freed of errors 
at any particular order in a by a perturbative caluclation with a finite number of diagrams. 
The only corrections which can occur at 0{a~^) are a correction to the mass, 

6mf ~ QcO''^ , AcO"^ (8) 

and a constant difference between lattice and continuum values for (0^)5 

~ a-' . (9) 

These are generated by one loop diagrams, and the latter is already present at zero coupling. 
Here and in what follows the 3 subscript will be implied, when there can be no confusion. 
At 0(a°) or 0(a° In a/i), the only corrections are of form 

Sm^r^g^^, gX, (10) 

and 

Ac. (11) 

The coefficients are determined by a two loop calculation and may contain ln(/ia). Knowing 
all of these corrections allows one to make a match which is correct in the strict a ^ limit. 
These counterterms are computed for SU(2) Higgs theory in [^, Q and extended to SU(A^) 
in p. 

At 0(a), the possible additive corrections to mf and ^'^^ 

6mf ~ ag^ , ag^X^ , ag^X^ , aX^ , 

and 

~ ctfi'c 5 ag^Xc, aXl . 
These arise from 3 loop diagrams. The other possible corrections are 

mf — ml ~ o-glml aXc-ml , 
(0V)«-(0V)c ~ agl{(f)'^(f))c, aXc{(j)'^(j))c, ami, 
Z<p-1 ~ agl, aXc , 

SXi ~ ag"^ , aglXc , aXl . 

These corrections are potentially generated by one loop diagrams. 

We should also check whether nonrenormalizable operators are induced; the terms we 
might possibly need to include in the tree action are of form 

{agl , aA,)(0V)' , aF^<f)^<j), a(t>^{D4f , a(0V)' • (19) 

But any diagram which generates one of these contains too many powers of g^ or A, which 
must be balanced by more powers of a. Even at one loop, none of these terms are generated 
below 0{a?). So the corrections listed are exhaustive at 0(a). 
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(12) 



(13) 



(14) 
(15) 
(16) 
(17) 
(18) 



It would be nice to compute all 0(a) corrections. But we see that two of these, additive 
corrections to mf and respectively, involve 3 loop computations, which are prohibitively 
difficult. Hence a complete analytic treatment of all 0{a) corrections appears unrealistic. 
Nevertheless, we accomplish something worthwhile by computing the remaining coefficients. 
This is because they influence and report the strength of the phase transition, while the 
additive corrections to mf and do not. Computing the others eliminates 0{a) errors 

in jumps in order parameters, the surface tension, the latent heat, the profile of the phase 
interface, the temperature range of metastabihty, etc. Only the MS value of mf^, at Tc and 
the absolute MS value of o^is phase are determined with less accuracy. Further, it 

may be much easier to determine the (uncalculated) coefficients for the 3-loop mf and 
corrections by Monte-Carlo methods if all other 0(a) errors are eliminated, and it may be 
easier to compute these than it would be to eliminate all 0(a) errors in all quantities by 
extrapolation. 

Finally we should point out that the 3-D values of m| and cire known with less 

accuracy in terms of the 4-D parameters (and hence physical measurables) as g"^, A3, and 
wave function normalizations. Writing 

gl = gl{fi)T (1 + cMT)gl + c,{ii/T)gl + ...), (20) 

a 1 loop calculation of the dimensional reduction step already establishes the value of Ci, 
and the relative uncertainty in g1 is 0(q;^). On the other hand, a 1 loop determination of 
'^L/fi'lc leaves an 0(1) error, and the 2 loop determination, which is presently state of the 
art, leaves an 0(a) error.0 Improving this by doing the 3 loop computation would be of 
comparable difficulty to the 3 loop calculation of the lattice continuum match. If we make 
the parametric estimate that gl^a ~ a, ie a ~ 1/T, then the 0(a) calculation of 5m\i is 
unjustified without the 3 loop dimensional reduction calculation. But the 0(a) calculation 
of the other quantities is parametrically justified. 

It is possible in principle to continue the lattice to continuum match to arbitrary order in 
a. At any finite order, the lattice action requires a finite number of terms and the matching 
requires the calculation of a finite number of diagrams. In practice the complexity rises very 
fast beyond 0(a). To remove all O(a^) errors in measurables describing the strength of the 
phase transition, one must expand the lattice action to contain terms which can be tuned 
to remove all tree level 0{a?) nonrenormalizable operators. This changes the Feynman rules 
of the perturbation theory. One will also need to compute corrections to the wave functions 
and couplings at two loops, which involves over 200 topologically distinct diagrams. We will 
content ourselves with 0(a) improvement here. 

3 Improvement in SU(A^) theory 

Our goal is, given a continuum SU(A^) plus real adjoint Higgs plus complex fundamental 
Higgs theory in 3 dimensions, to write down the lattice theory which is equivalent up to 0(a) 
additive errors in masses and expectation values of operator insertions and 0{a?) errors in 

•^This is for the integration over the "superheavy" scale, ie nonzero Matsubara frequencies. If one also 
integrates out zero Matsubara frequency modes which take on 0{gT) masses, the "heavy" scale, then a in 
this argument becomes a^^^. 
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wave functions, couplings, multiplicative corrections to masses and operator insertions, and 
nonrenormalizable operators. 

We denote the adjoint Higgs field as $ = ^"■T'^ and the fundamental Higgs field as (p. Our 
group conventions are standard, ie TtT"'T'' = {l/2)6ab- It is to be understood that Di acting 
on $ is the adjoint covariant derivative and Di acting on (p is the fundamental covariant 
derivative. The Lagrangian density of the continuum theory we consider is 



,ont = -TrF,,F,^+Tr(A<fA<f) + (A0)^(A 
+m|Tr$^ + m 



Ai(Tr$^)^ + A2W + A(0V)' + 



(21) 



Note that there are two possible adjoint field self-interactions and two possible interactions 
between fundamental and adjoint Higgs fields. In SU(2) and in SU(3), Tr$^ = (l/2)(Tr$^)^. 
In SU(2) this is because dabc = 0, and in SU(3) it is because of the relation, special to SU(3), 
that 



1 



{SabScd + SacSbd + SadSbc) for SU(3) . 



(22) 



Hence, in these theories the two adjoint Higgs field self- interact ions are dependent, and only 
the combination (Ai + A2/2) is important. Also note that the /i2 term can be dropped for 
SU(2), because dabc = in that case. 

A choice for the Lagrangian density for the lattice theory which is general enough to 
allow 0(a) corrections is 



latt (a;) 



2Z 

(X 



2Zi 



^[Tr$2(x) - Tr$(a;)[/,(x)$(x + i)U}{x)] + 

i 

+(A + (5A)Z|(0t0)2 ^ ^ <5Ai)Z|(Tr<l>2)2 + (A2 + 5A2)Z|Tr<l'4 + 
+ {hi + 5/ii)Z^Z$0VTr$' + ^^^^^^Z^Z$rf,6,<l>"<l>V^T^0 + 



+ 



Z$Tr$^ 



(23) 



The couplings and masses here are the continuum theory values, and the coefficients Zg, Zq>, 
Z^, 6Xi, Smf, and Z^ perform the 0(a) corrections. Z^ is a 2 x 2 matrix, which differs at 
0{a) from the identity matrix. 

Further, to convert measured values of operator insertions to the equivalent continuum 
values, we need an addative and a multiplicative renormalization of the operator insertions. 



(Tr$2) 



continuum 



5(0t0) 

(5(Tr$2) 



+ Z 



OP 



(Tr$2) 



(24) 



latt, measured 



Here Zqp is also a matrix, with 0{a) off diagonal elements which account for mixing between 
the operator insertions. All notations are chosen to follow [|19]; to convert to the notation 



7 



of ||T3[, change the signs on the counterterms in the scalar effective potential, introduce 
Xla = 4:Xi/g'^, rescale the scalar fields by a factor of g'^aT/4 and the by a factor of a^, 
and write A/ g'^a = Pl- Also, Zg was called Z^^ there. 

Neither Z^j), nor Zqp are separately gauge invariant or physical; one may for instance 
rescale the $ field by a factor of r/ = 1 + 0(a) and change Z$ by rj^'^ and the right column 
of Zqp by 7]^"^ without changing anything. It is the matrix Zqp x Diag[Z^^ Z^'] which 
is gauge invariant. This combination equals Z^, a statement which will remain true at all 
orders. 

Similarly, while Zg is gauge invariant, if one tries to relate the 3-D theory gauge field 
A to the lattice gauge field A defined by Ui{x) = exp{—iagT°- A'} {x + one will find 

the relation contains a non-gauge invariant 0(a) correction, for the same reason that the 
4-D continuum gauge field wave function renormalization is not gauge invariant, while the 
renormalization of the coupling constant is. This is only important if one is interested in 
measuring A field correlators in some gauge, and we will not compute the A field wave 
function correction here. 

To determine all of the required corrections, one must calculate all 1 loop self-energy 
insertions for all fields to O(p^), all 1 loop 0{p) corrections to one scalar-gauge vertex, all 
1 loop scalar 4-point corrections at zero external momentum, all 1 loop corrections to the 
insertion of a 0V or Tr$^ operator on a zero momentum scalar line, all 1 and 2 loop scalar 
self-energy insertions at zero momentum, and all 1 and 2 loop vacuum diagrams with 0^0 or 
Tr$^ insertions. The topologically distinct diagrams required in Landau gauge are presented 
in Figure |I]. In a general gauge, additional diagrams are needed. The Feynman rules are 
spread between ||, |1^ 



All of the diagrams have been computed previously: (a) at 0(p°) in [^, (6) in [^, (c) in 
0],"^ and (a) at O(p^), (d) in [0. Our task here is purely combinatoric. We must include 2 
types of scalars and compute the appropriate SU(A^) group factors. All the group theoretic 
identities which are needed appear in Appendix A of [|1^]. The converse is also true. 

Rather than present the full combinatorial details of the calculation, we will skip to 
the results. Two numerical constants appear, and their values are S = 3.175911536, ^ = 
0.152859325. Also, when comparing the dimensionality of the right and lefthand sides, 
remember that the coupling constants g"^, Aj are in 3-D notation and correspond to T times 
the dimensionless 4-D values (at some renormalization point). The coupling renormalizations 
are| 

V- - $(|^(2--3)..(f -^).4il4lWl^-|IV (25) 



6\ 



a 



'N^ + N^-4:N + 2 




An \ AN^ 



*The computation there is not organized in terms of two loop scalar self-energy insertions, but it is 
equivalent. 

^We choose to report Z^^ — 1 because at two loops this quantity is only changed by two loop self- 
energy diagrams and vertex corrections, ie it absorbs interated insertions of 1 loop self-energy diagrams on 
propagators. But at 0{a) accuracy we could equally have written Eq. ( p5| ) as a correction to 1 — Zg. 
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4AA2 — fi 

+^^A^e + ^hl^ - ^hl^ - jm + S)/A,) , (27) 

/V2 _ Q 1 \ 

^2 = ^[Ng'^ + 2^;^Xl^ + 12X,X,^+-hl^-j{m + j:)g'X2j, (28) 

7V^ — 4 

= ^(2/e+(A^' + l)/^iAie + (2Ar + 2)/iiAe + 2/i^e+ ^ /^2^ + 

Sh2 = ^ (^Ng'^ + h2 (^2X + 2Ai + ^^^^Aa + 4h + ^2iV^^ ^^) ^ ~ 

S Af2 _ 1 \ 

:i8e + S)//^2 . (30) 



12Ar 

For the theory without the adjoint Higgs field, drop the term with square brackets in Eq. 
(p5|), set Ai = A2 = /ii = /i2 = in the equation for 6X, and ignore the equations for 6X1 
etc. Putting = 2, we recover the results of [|l^, though the notation there is somewhat 
different. For the theory without the fundamental Higgs, drop the term in curly brackets in 
Eq. (p5D , set X = hi = h2 = 0, and ignore the equations for their corrections. 

The fundamental scalar contribution to Eq. ( P3| ) is wrong in by a factor of 1/2, both 
for SU(2) and U(l). This error is not numerically important; even in SU(2) it leads to a 3% 
error in Z^^ — 1. Appendix C of that paper, which treats the case with both Higgs fields 
for N = 2, has several typographical omissions of factors of S/47r and ^/47r, and also the 
g^ contributions to {6X1 + 6X2/ 2) and to 6hi are both off by a factor of 2 there, due to a 
combinatorial error. These are the only combinatorial mistakes in that paper.Q 

Before continuing, we mention two checks on the calculation so far. First, both for N = 2 
and = 3, the combination (Ai + A2/2) is important, not Ai and A2 separately. Checking, 
one finds that {6X1 + 6X2/2) is fixed when one varies Ai keeping (Ai + A2/2) fixed, both for 
A^ = 2 or A^ = 3. Also, /12 should have no infiuence when N = 2. Indeed, its contribution to 
the rescaling of other couplings vanishes for N = 2. 

Now let us present Z^, and Zqp- Since they are not separately gauge invariant, 
we will present them in Landau gauge. In this gauge, at one loop, the corrections to wave 
functions all go as g'^ and the corrections to Zqp depend on scalar self-couplings. This 
statement does not persist at 2 loops. The results are 



An 



{2N + 2)A Nhi 

(N' + l)Ai + ^A2 



(33) 



^However, early preprint versions had further mistakes which were corrected before pubhcation. 
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Here 1 is the 2x2 identity matrix. Varying the gauge changes the wave function corrections 
and the on diagonal elements of Zqp, but not Z^- The off diagonal elements of Zqp are also 
gauge invariant at 0(a), but this will not persist at higher loop orders. 

It remains to present Sm"^, 5m|, 5{(j)^(j)), and 5(Tr$^). The one loop corrections are [jl9 

1 



5m0(one loop) 
5m\ (one loop) 
5{(j)^ (j)){one loop) 
(5(Tr$^)(one loop) 



Ana 
S 
47ra 

-A 

Ar2 



y + 2(iV+l)A + 



-h 



2Ng^ + Nhi + (A^2 + 1)^^ + 



2N^ 



N 



-A, 



Ana 
- 1 



An J ' 
47r 



47ra 



(34) 
(35) 
(36) 
(37) 



We have included the 0(a), dependent correction to the operator insertion expectation 
values. For the case where there is one Higgs field but not the other, set hi = 0. Also note 
that our hi differs by a factor of 2 in normalization from the parameter h^ appearing in [^. 

The two loop corrections have been computed in |T^, for the case where there is an 
adjoint scalar, or a fundamental scalar, but not both. For our case, we must add new terms 
when there are two types of scalars, and contributions from inserting 0(a) counterterms into 
1 loop diagrams. 



The two loop contributions listed in |]I9 



+ 



5(Tr$2) + 



167r2 
+2\g 



9 TTTo + 2A^7' 



,4iV^ -iV + 3 



N 



4Ar2 



4Ar2 



AN^ 



aX' 



1 



4 



-T? + 



A^- 1 
3Ar2 



8 o 
-vrS + iV^ + 



+ 1 - 4A^(A^ + l)p - 2(3A^2 - 1)5 + 2N^{2ki - k^) 



(38) 



167r2 



2iV(72 (Ar2 + i)Ai + 



2A^2 



4(2A^2_3) jy4 

-A1A2 



6Ar2 _ 



-A2 - 2(Ar2 + l)Xl 



N 



ATS 



-A 



. a/i / 



+2Ar/ — -5 X (Ar2 + l)Ai + 



^ A 



2N^-3 
N 



A2 + 



+ 



3N^ -8 
6A^2 



-TT 



S-4(5 + p)+2ki-K4 



167r2 
,2 



\ 6 , 
In— + C + — 
a/i 4 



5 



IGvr' 



-A^(A^2 



',6 ^ 

In — + C + — - 5 
y a/x 4 y 



(39) 
(40) 
(41) 
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where the newly introduced constants are ( = .08849, S = 1.942130, p = —.313964, ki = 
.958382, K4 = 1.204295. Here and in the following we will write + = to show that these 
contributions are to be added to those listed previously. 

In the case where both Higgs fields are present, there are added contributions to the mass 
counterterms of 



6ml + 



6ml 



+ 




hi + Nhig"^ 



1 4 
4^ 



InA + C^ 
, a/i - 



+ 



P9 



(42) 



hi + (iV2 
- l)h^g^ - 2iVp^?4 




(43) 



At = 2 these match the results of 0]. Do not include these in the theory with only a 
fundamental or an adjoint Higgs, but not both. 

In addition, there are contributions at 0{oP) arising from 0{a) conterterm insertions into 
diagrams which give 0{l/a) divergences. These lead to the following corrections: 



5ml + 



6ml 



+ 



6{Tt^^) 



Ana 

S 
Ana 

y _ 



N 



1) + 2{N + l)6\ + 



1 



6hi 



2Ng\Zg 

5(0t0) 
5(Tr<l>2) 



1) + N6hi + (A^2 + i-^^Xi + 
1 loop) . 



2A^2_3 
N 



(44) 
(45) 
(46) 



The 0{a) counterterms appearing here were given earlier, and as usual the term in the 
parenthesis in Eq. ( ^BD is to be interpreted as a 2 x 2 matrix. It should be clear how to pare 
these equations down to the case where there is only one type of Higgs field. 

This completes the calculation of all counterterms needed for the 0{a) improvement of 
the theory, except for the additive 0(a) corrections to the masses and operator insertions 
which arise at 3 loops, which we will not calculate, as advertized. 

What should you do, though, if you have already taken data without including these 
corrections in the Lagrangian? That is, what if you have data for a theory which naively 
has a lattice spacing a„, no wave function or g renormalization, and naive couplings A„, etc? 
Then you should figure out what theory you were "really" looking at, and reinterpret the 
results accordingly. By demanding that the "naive" action 

«n E I ^ E Tr(l - P^Ax)) + 4 Y.[^i{x)^n{x) - 0|,(x)f/,(x)0„(a; + l)] + 

X [9 (^n i<cj i 



+ - EtTr-^"' (a;) - Tr$„(x)f/,(x)<l>„(x + t)u!{x)] + 



a,^ 



+ Ai„(Tr$^)^ + A2.Tr$ 
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(47) 



be equal to a J2x ^lattix) of Eq. (pSf), one finds that the naive values labeled with the n 
subscript are related to the 0(a) corrected values by 



g 

(ZgZ, 



1^1/2, 



A 
Ai 
A2 
hi 
h2 
5ml 

6 ml 

ml 
ml 



ZnXlr 



ZnX2n 



5\i, 
Zghin — 5hi , 

Zgh2n — 5h2 , 

Zg5ml „^ , 
ZgSml^^ , 



Z'gZ-' 



m" 



(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
(54) 
(55) 
(56) 
(57) 

(58) 



Note that, up to 0(a) corrections, 5ml computed using a„ and the naive scalar self-couplings, 
and without including Eq. (i3), differs from 5ml computed using the improved values and 
including Eq. (|4^) , by precisely a factor of just what is required above.[] This just says 



I9| is the right one if you use the "naive" 



that the mass renormalization computed in 
couplings and wave functions. 

Next we should relate the 0(a) corrected value of the continuum operator insertion 
expectation values, through the 0(a) corrected lattice values, to the uncorrected lattice 
values. We have 



(Tr$2) 



continuum 



5(Tr$2) 



+ z 



OP 



ZgZ^ 



(Tr$2) 

+ (5(0V))n 

+ (5(Tr$2)), 



latt, measured 



(59) 



Again, using a„ in Eqs. ( |36| , |37D and skipping Eq. ( |46D leads to a value for the naive 
counterterm which matches the corrected counterterm in the above. So if one has "naive" 
measured values for condensates, with the additive counterterms from ^, |T9[ already 
taken out, then one applies the 0(a) corrections by multiplying by ZgZ^, where Zm is to be 
understood as a 2 x 2 matrix, as usual. 



^To see this, note that one factor of Zg conies from using a rather than a„ in Eqs. ( |34| , |35| ) and the other 



comes because Eqs. 



5|) multiply the gauge part by Zg, while the scalar part goes as A + ^ 



ZgXn 
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So to apply 0(a) corrections to data taken using a "naive" lattice Lagrangian, interpret 
the true lattice spacing as Z~^a„, the true scalar self-couplings as ZgXn — dX, the true mass as 

= Z^^ZgiTin, and the true value of order parameters as Z^Zg times the naive values. Note 
that Zg < 1 and {Zg — 1)A — 6X<0, so jumps in order parameters are generally smaller than 
the "naive" values and the true scalar self-couplings are generally smaller than the "naive" 
values. Using the "naive" values tends to over-report the strength of the phase transition. 
Also remember that the mass correction should only be used on the difference in 3-D mass 
parameter between two trials with the same lattice spacing and scalar couplings, and the 
order parameter correction should only be used on the difference in the order parameter 
between phases or between trials with different 3-D masses (and here one should also be 
careful to include the mass dependent term in Eqs. ( |36| , p7D ). This is because we have not 
computed the 0(a) additive corrections to the mass and order parameter counterterms. 

To correct the latent heat and the surface tension, note that they are energies per unit 
volume and area, respectively; the 0(a) corrections come about from the relation between 
a and a„. Similarly, correlation lengths are corrected by converting from lattice units to 
physical units using a rather than a„. If one is interested in operator insertions not discussed 
here, it will be necessary to extend this work by computing the renormalizations of those 
insertions. 

We will present numerical values of renormalization constants in some particular theories, 
and give step by step instructions for building 0(a) improved Lagrangians, in Appendix B. 

One final comment is in order. As we have written them, the corrections Zm, Zg, 6X etc 
which are needed for this 0(a) improvement are to be computed with the 0(a) improved 
parameters A, a, etc, not with A„, a„, etc. If one only knows the "naive" values, this might 
require an iterative or "bootstrap" type calculation. However, the difference between using 
the improved and unimproved parameters in the calculation is formally O(a^), so one should 
not worry too much which is used-except for computing the one loop 0(l/a) corrections to 
masses and order parameters, as already discussed. 



4 What to do in SU(iV)xSU(M) 

Sometimes it is interesting to treat a case where the gauge group is not simple. Here we will 
restrict our attention to SU(iV) xSU(M), where each scalar field is a singlet under one or the 
other group, and we intend to be illustrative but not exhaustive, so we will consider the case 
where an SU(A^) group has both an adjoint and a fundamental scalar (it is easy to delete one 
of them) and there is one more scalar which transforms trivially under SU(A^) but perhaps 
nontrivially under the other group, and has self-interactions. Calling the added scalar s, and 
writing to mean 2s^s if it is fundamental and 2Trs^ if it is adjoint, the new interactions 
between the s field and the scalar fields and $ which can appear in the Lagrangian are 

^s^^^+^s^Tr^\ (60) 

If hsi = = hs2 then the theory decouples into two noninteracting parts, the SU(A^) part 
and the part which contains the scalar s. We assume that one can compute the corrections 
in this limit by using the last section, and we will compute the new corrections which must 
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be added at nonzero hgi or hs2- We will compute the corrections to the couplings of the 
SU(A^) sector and hgi, with the minimum information about the scalar s. To compute 
corrections to the sector containing the s field, just think of it as the SU(A^) sector and use 
these same results. 

We need only a httle information about the scalar s and its (gauge and self) couphngs 
to compute the corrections to the SU(A^) sector. What we need to know is the number of 
degrees of freedom s contains, A^dof (in SU(M) this is 2M for a complex fundamental field 
and — 1 for a real adjoint field), and its one loop mass renormalization at hgi = hs2 = 0, 

6ntt(one loop , zero hsi, hs2) = --, — (Ci + C2) , (61) 

47ra 

where Ci is the contribTition from all scalar couplings and C2 is the contribution from any 
gauge couplings it participates in. They have the same units as a coupling constant. 
The one loop mass corrections due to the new particle are 

6m, -4^^/^.!, (62) 

X 2 , ^ ^DOF, 

6m^ + = — hs2 ■ (63) 

47ra z 

there are no new one loop additive operator insertion corrections. 

The one loop 0(a) corrections to scalar self-couplings are increased by 

(64) 

(65) 
(66) 
(67) 

(68) 

^—h^hs2 + 2hi^ i - 



sx 


+ = 


Ndof ^2 
Att 4 ^^1' 


5Xi 


+ = 


{a A^DOF ^2 
Att 4 


5X2 


+ = 


0, 


5hi 


+ = 


Ndof , , 
47r 2 


Sh2 


+ = 


0, 


Shsi 




(c,h,i + 2{N+l)Xh,i + 



^^.^ + §)(18e + E)/..,j, (69) 
Shs2 = £| {C,hs2 + {N^ + l)Ai/i,2 + ^^^~^ A2/i.2 + Nhhsi + 2hl^^ e - 

f/ + §)(18e + S)/.32}. (70) 

There are no corrections to the gauge field renormalization Zg at this order. Such corrections 
actually start at O(a^). 

Moreover, we should now think of Zqp as a 3 x 3 matrix, which at hsi = = hs2 is 
nonzero in the upper 2x2 block, where it equals Zqp of the last chapter, and the lower 
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right component, where it is Zqp for s. There are no new contributions in these bfocks. The 
terms we should compute are 



Ana 



same same Nhgi 
same same — 7^ — h 

— 5 — ''-si — 9 — '''s2 



2 ■-52 

same 



(71) 



where "same" means the same value as at hgi = = hs2- At one loop there are no corrections 
to wave function renormalizations due to s. These will first appear at O(a^). 

What remains is to compute two loop mass and operator insertion additive corrections. 
The only new correction to the operator insertion conterterm is that the new value of Zqp, 
Eq. (ItTD , should be used in Eq. (^). The new corrections to the mass counterterm, aside 
from using the new values for 6X etc. in Eqs. p^j45|), are 



Sml + 



AT, 



DOF 



2 

A^DOF 



dhsi + 



Area 
S 
4:7ia 



In — + C 1 + C2 
a/i / 




C, —-6 



,(72) 
.(73) 



If there are more than one scalar which are SU(iV) singlets, the corrections are just the 
sum of the corrections from each particle, except for cross terms in the renormalization of 
the hg type interactions, which can be figured out from Eqs. ( |69| , |70D . Things become more 
complicated at higher loop order, however. If there are 3 or more gauge groups, each with 
scalars, then the combinatorics of the one loop scalar coupling corrections get messy and 
what we show here may not be sufficiently general. But we know of no physically interesting 
theories in this class. Also, if there are scalars which transform nontrivially under more than 
one gauge group, the situation becomes considerably more complicated. For instance, the 
one loop self-coupling and two loop mass counterterms then contain gfg2 type terms. The 
particularly simple case of SU(2)xU(l) is treated in p3| , |T9[. We will not attempt anything 
more general here. 



5 Conclusion 

We have discussed why it is possible to write 0(a) corrected actions for 3 dimensional lattice 
Yang- Mills Higgs theories, and we have extended the previous results for SU(2)xU(l) to 
SU(A^) with a fundamental and/or an adjoint Higgs field. We have also shown how the 
corrections can often be applied to uncorrected data "after the fact." Generally the result is 
that the actual value of 3-D scalar self-coupling used is lower than the naive coefficient put 
into the lattice calculation, and actual jumps in the order parameters are smaller than the 
uncorrected data imply. 

The limitation of the work here is that we have not been able to compute the 0{a) 
additive corrections to the Higgs field masses or the 0^ operator insertions. Computing 
these is possible in principle but requires doing a lot of very nasty 3 loop diagrams. We 
have argued that for many purposes these are the least important 0(a) errors, and that the 
knowledge of the relation between the 3 dimensional and physical, thermal 4 dimensional 
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theories is in any case weaker here. The corrections for the other coefficients are definitely 
valuable for numerical efforts to study hot gauge theories nonperturbatively, allowing better 
precision with less numerical effort. 

A final issue we should address is: how small should a be made so that we don't have 
to worry about the remaining O(a^) errors? Of course, the answer depends on what one is 
doing. For instance, if one is studying the strength of a phase transition where it is strong, 
say at very small A or Ai, then the natural scale of physics involved becomes shorter than 
1/ and nonrenormalizable operators may be a problem. One can estimate the importance 
of these operators by including them in a perturbative calculation of the effective potential, 
as was done in |]13|. The conclusion is that nonrenormalizable derivative terms are less 
important than the 0{a) corrections 5 A to the scalar self-couplings, roughly by a factor of 
ga\J (p'^cp. If neglecting the 0(a) scalar self-coupling corrections would make a significant 

difference and ga^J^Up is not very small, then a is not small enough. 

Away from the case where the phase transition is strong, the two loop, O(a^) corrections 
to the Lagrangian parameters may also be important. We expect that the largest of these 
will be "tadpole" type corrections to the gauge field normalization, higher loop analogs of 
the first term in Eq. (pS]). This term arises from a nonrenormalizable operator in the action 
of form TrF-^-. It turns out to dominate Z^'^ — 1, which in turn is the largest of the 0{a) 
corrections. We might expect that "tadpole" type corrections, from multiple appearances of 
this term and from the term of form TrFj^j, give the most important contributions to — 1 
at two loops, and that this is again the largest 0{a^) correction. (This would follow a fairly 
general pattern of behavior for lattice gauge theories.) Fortunately, it is not too hard to 
compute this limited class of 2 loop diagrams; we do so in Appendix A. The conclusion is 
that the two loop correction to Z^^ — 1 is about the same size as the square of the one 
loop correction. So if the one loop correction is 10%, then the two loop correction is 1%. If 
nonrenormalizable operators are not a problem, then the one loop improvement is likely to 
be sufficient as long as its coefficients make corrections in the 10% range. 

It may not always be possible to reduce O(a^) errors to the level desired, and sometimes 
it is also important to know mass parameters with high precision. In these cases it may still 
be necessary to perform an extrapolation over data at several values of a. We expect that the 
0{a) corrections will still be useful in this case, because they should make the extrapolation 
much better behaved. 
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A Two loop tadpole type diagrams 

In this appendix we compute two loop "tadpole" corrections to the gauge field self-energy. By 
"tadpole" corrections we mean corrections emerging from a certain restricted subset of the 
interaction terms generated by the gauge field term in the action. Namely, we approximate 
= exp{iga'^T"'F°'-), where Fij is linear in the gauge field, ie in Fourier space it is 
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F^j{k) = kiA'j{k) — kjA°;{k). This amounts to neglecting appearances of fabc in favor of Sab 
and dabc- 

The terms in the action which emerge from this approximation are 



a 



These give the gauge field propagator and a tower of interaction terms. The term with 
2n powers of Fij gives a vertex which attaches to 2n lines, with momenta we will label as 
Pi, p2, . . ., polarization indicies we label as ii, Z2, . . ., and group indicies we label as a, b, . . .. 
The Feynman rule for the vertex is that it contributes 



jk 

J_ (^Tr(T''T* . . .) + all permutations of a, 6, . . .) . (75) 



(2n)! 

Though the nonrenormalizable terms here have several powers of a, they lead to diagrams 
with strong power law UV divergences and are important to the radiative corrections. The 
first of these is responsible, through diagram (a) in Figure |^, for the first term in Eq. (p5|). 

The first natural question is, is this approximation any good? To answer that, compare 
the numerical value of the first term in Eq. ( |25D to the sum of the others. Even with both 
scalars present, in SU(2) the first term is larger than the others combined by a factor of 
3. At large the domination is by a factor of 6. Without scalars the domination is twice 
as large. So it appears that "tadpole" type terms do give the dominant contributions. We 
should also comment that the contributions from these diagrams are gauge invariant and 
transverse diagram by diagram. 

Now let us compute the two loop "tadpole" type terms. There are insertions of the 0(a) 
counterterm in the propagator and the vertex of diagram (a), but these cancel. There are 
also diagrams (6), (c), and (rf), which give contributions to Z^^ — 1 equal to 

c/Vtt^ (N^ - 6N^ + 10\ /2 X 0.3015814\ 

^ [ em ) [ 3 ) ' 

- ^ (^) Q ^ 

We have written each contribution as a common "loop counting" type term, a group theoretic 
term (with powers of g'^/4: taken out into the leading expression), and the product of the 
symmetry factor and momentum integration. Diagram (c) gives precisely the square of the 
one loop contribution. For every value of A^, the sum of these 2 loop diagrams is less than 
but on order the square of the one loop correction. 
The numerical constant, 0.3015814, is the integral 



I 



d^kdHd^ 




{27T)'^6\k + I + m) , (79) 
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where ki = 2sin(fcj/2) and fc^ = J2i kf- This integral, hke all the integrals arising from the 
tadpole terms, is completely infrared safe and has no continuum anolog. 



B Some particular examples 

In this appendix we will plug in numerical values and give step by step instructions for 
building 0{a) improved Lagrangians for the specific theories of most interest, namely SU(2) 
plus fundamental Higgs, SU(3) and SU(5) plus adjoint Higgs, and SU(3) xSU(2), each with 
a fundamental Higgs. The latter is a model for studying the phase transition in the minimal 
super symmetric standard model with a light stop squark, which makes including the SU(3) 
color sector necessary. 

In each case the first step is to write down the desired values for g^a/A, Xi = Xi/g^, and 
yi = mf/g'^, which one wants the simulation to correspond to. Then one computes a number 
of counterterms and uses them to find the "naive" values 5'^a„/4, Xi^n, yi,n corresponding to 
the desired improved values. One constructs the lattice Lagrangian using these, computing 
two loop mass squared and operator counterterms and 5(0^) exactly as one would in 
ignorance of the 0(a) improvement scheme, using expressions in |], [19|. Then one applies 
0(a) multiplicative corrections to "naive" measured values of ((0^) + 5(0^))naive and the 
surface tension (Xnaive- 

B.l SU(2) + fundamental Higgs 

We begin with SU(2) plus fundamental Higgs theory. Writing x = Xjg^ and y = rn? /g'^, we 
choose desired values for g'^a/A, x, and y, and then compute the counterterms 



2 

Z;'-l = 0.6674^, (80) 

5x = (0.01825 - 0.4717a; + 0.5839x2) 9^ ^ ^g^^, 

Zm-l = (-0.2358 + 0.2919a;) ^ . (82) 



Now compute the "naive" values 



g ttrr = Zgg a, (83) 
= Z;\x + 6x), (84) 



X. 



Vn = Zg'^Zmy, (85) 

and construct a lattice lagrangian which, according to naive lattice to continuum relations, 
should model the continuum theory with these "naive" values. Compute the counterterms 
5yn and {6{(f)^(f)))n using the two loop expressions in 0, and using a„, x„. Do not include 
the corrections, Eqs. (H3,H5|H3). (Alternately, compute them using g^a and x, and include 



Eqs. ( ^^^6|) , but then multiply by Z '^Zm- The two approaches differ at 0{g^a) and it 



is not clear without the 3 loop calculation which will result in smaller 0(a) additive errors.) 



18 



This theory will in fact represent the continuum theory with lattice spacing a, scalar 
self-coupling A = xg^, and Higgs mass = yg'^, up to 0{a?) errors (and an 0(a) additive 
shift to y). The 0{a) corrected value for the order parameter is 



-'^0) corrected 



+ (5(0V))n) 

and the surface tension is 

C"corrected -^(^^naive • l^"^) 

Note that in 3-D units the surface tension is just an inverse area. To convert to 4-D units 
one multiplies by T. Physical lengths such as correlation lengths should be corrected by 
multiplying by Z^^. 

B.2 SU(3) + adjoint Higgs 

Suppose we want an improved lattice action for SU(3) plus adjoint Higgs theory at a desired 
value of 5'^a/4, x = (Ai + A2/2)/(7^, and y = m\/g^. We compute 

2 

Z-^-l = 1.3299^, 

9"a 



6x = (0.21895 - 1.8867X + 0.7785a;2) ^ , (89) 



4 

2 

Zm-l = (-0.9434 + 0.4866a;) ^ . (90) 

Note that the pure gauge (no powers of x) corrections are much larger here than in SU(2) 
plus fundamental Higgs theory, particularly for the self-coupling. This means we need a 
smaller value of g'^a/4 to fight down lattice artifacts, especially if we are exploring small x. 

From this point the steps are the same as in SU(2); write down a lattice action which 
according to "naive" lattice to continuum matching corresponds to the continuum theory at 
the "naive" values 

g'^an = Agg^a, (91) 
Xn = Zg\x + 6x), (92) 

Vn = Zg'^Zmy, (93) 

and compute the counterterms for the mass squared and operator insertion from the equa- 
tions in |]T^. Take data, and correct the Tr$^ order parameter by multiplying by ZgZ^ and 
the surface tension by multiplying by Z^. 

B.3 SU(5) + adjoint Higgs 

This case differs from SU(3) plus adjoint Higgs theory only in that there are now two scalar 
self- couplings Xi = \i/g^ and a:2 = ^2/9"^, and that there are different numerical values for 
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the counterterms, 



Z,^-l = 2.454^, (94) 



= (0.1460 - 3.145x1 + 1.557x? + 0.915x1X2 + 0.163x^) (95) 
5x2 = (0.2433 - 3.145x2 + 0.584x1X2 + 0.311x^) (96) 

2 

Z.^-l = (-1.572 + 1.265x1 + 0.457x2)^. (97) 

From here on the procedure follows the SU(3) plus adjoint Higgs case. Note that g^a/A must 
be even smaller for SU(5) than for SU(3), which is not too surprising. 

B.4 SU(3)xSU(2) + fundamental Higgs in each 

This theory could be phenomenologically interesting if the left handed stop squark is light 
enough that a perturbative treatment becomes unreliable. It holds a potentially rich phe- 
nomenology, allowing a double phase transition and an exotic color broken symmetric phase 
|T^. It also nicely illustrates what to do when there are more than one gauge group, and 



it presents a heirarchy problem which makes its numerical study without 0{a) corrections 
almost impossible, and its study even with them rather tenuous, for the physical value of 
9s/9l- 

We will use the notation of Laine and Rajantie. The continuum Lagrangian is 

£ = ^TrGj + iTr4 + (Af/)^(At/) + (Ai/)^(Ai/) + 

+mlU^U + mjjH^H + Xu{U^U)^ + XniH^Hf + -fU^UH^H , (98) 

where U, the squark field, is in the fundamental representation of SU(3) (the field G, with 
coupling constant gs) and H, the light Higgs field, is in the fundamental representation of 
SU(2) (the field F, with coupling constant g^u)- 

The theory is characterized by one dimensionful number g"^ and 6 dimensionless param- 
eters, xi = Xu/gs, X2 = Xn/gl, = 7/^s> yi = ^u/at, 1/2 = mjj/g^, and z = gl/gl Note 
that we express all of the masses and couplings in terms of g^. To put the theory on the 
lattice, we first choose a desired value for each dimensionless number and a value for the 
lattice spacing in terms of the one dimensionful number, that is we must choose g1a/A. The 
0(a) improvement is carried out by first computing the counterterms 





- 1 = 


1.2619^, 


(99) 


z-J 


- 1 = 


(0.6674)z^ , 


(100) 




5xi = 


(0.03514 - 0.8386x1 + 0.6812x^ + 0.0487x^) ^ , 


(101) 




5X2 = 


(0.01825^2 - 0.4717^X2 + 0.5839x^ + 0.0730x^) ^ , 


(102) 
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^-0.4193x3 - 0.2358ZX3 + 0.3893xiX3 + 0.2919x3X3) 



9la 



-0.4193 + 0.3893x1 
0.0973x3 



0.1460x3 
-0.2358Z + 0.2919x2 



(103) 
(104) 



where we choose the convention that U'^U goes in the upper and H^H in the lower row in 
expressions where Zm acts on columns. 

Now we implement the 0(a) improvement by defining a set of "naive" paramters, 



yi,n 

y2,n 



Zgtt , 

Zs{xi + 6x1) , 

Zs{x2 + 5X2) , 
Zs(x3 + 5X3) , 

2 



z:^z 



yi 
2/2 



(105) 

(106) 

(107) 
(108) 
(109) 

(110) 



We then use a "naive" continuum to lattice relation to convert these into a lattice action, 
computing the mass and operator insertion conterterms from the expressions in [0. Finally, 



improved values for jumps in order parameters are computed by multiplying the unimproved 
values by the matrix Z^Zm, the surface tension is improved by multiplying by Z"^, and 
correlation lengths are rescaled by Z~^. 

One point is in order, though. For T ~ lOOGeV, the values of the coupling constants 
are around| g^jAn ~ 0.087T and gl/An ~ 0.032T. Hence, z ~ .368 ~ 1/3, so the "natural 
scale" of the strong and weak sectors differ by a factor of 3. To contain O(a^) errors in the 
strong sector, we need gla/A to be fairly small-for instance, we probably want Z'^^ — 1 to 
be on order 0.1. However, this forces a value of (5^^0/4) < 1/30, which is a very fine lattice. 
Yet we must make the lattice volume large enough to contain the rather long correlation 
lengths of the SU(2) sector, which puts some strain on numerical practicality. Hence there 
may realistically be something of a heirarchy problem involved in studying this system at 
the physical values for the couplings. 
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Figure 1: The topologically distinct diagrams needed for the Landau gauge calculation. 
Wavy lines are gauge propagators, solid lines are scalars, and dotted lines are ghosts. A 
cross is a self energy correction from the measure and a blot is a (fP' operator insertion. 
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(c) (d) 

Figure 2: Tadpole type diagrams which correct the gauge field self-energy at one and two 
loops. 
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